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Using event-by-event viscous fluid dynamics to evolve fluctuating initial density profiles from the 
Monte-Carlo Glauber model for U+U collisions, we report a “knee”-like structure in the elliptic flow 
as a function of collision centrality, located around the 0.5% most central collisions as measured by 
the final charged multiplicity. This knee is due to the preferential selection of tip-on-tip collision 
geometries by a high-multiplicity trigger. Such a knee structure is not seen in the STAR data. This 
rules out the two-component MC-Glauber model for initial energy and entropy production. Hence 
an enrichment of tip-tip configurations by triggering solely on high-multiplicity in the U+U collisions 
does not work. On the other hand, by using the Zero Degree Calorimeters (ZDCs) coupled with 
event-shape engineering such a selection is possible. We identify the selection purity of body-body 
and tip-tip events in full-overlap U+U collisions. By additionally constraining the asymmetry of the 
ZDC signals we can further increase the probability of selecting tip-tip events in U+U collisions. 


I. INTRODUCTION 

High energy collisions between heavy ions are used to 
probe emergent phenomena in Quantum Chromodynam¬ 
ics (QCD), the theory of the strong interaction. One 
feature of QCD is the transition from hadronic matter to 
a color-deconfined quark-gluon plasma (QGP) [1-3] as 
the temperature is increased. This transition can occur 
in heavy-ion collisions with sufficient energy for the sys¬ 
tem to melt into a hot dense fireball of asymptotically 
free quarks and gluons. 

Relativistic hydrodynamic models have been success¬ 
ful in describing the dynamical evolution of QGP [4]. 
Motivated as a testing ground for these models, a U+U 
collisions program was recommended in order to study 
the unique collision geometry resulting from the prolate 
deformation of the uranium nucleus [5-10]. Such a pro¬ 
gram was carried out in 2012 at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Lab [12]. 

To understand the attraction of uranium, consider that 
the initial temperature distribution of each QGP droplet 
is controlled by two main factors: deterministic collision 
geometry (i.e. the shape of the overlap region between 
two nuclei), and quantum mechanical fluctuations in the 
nucleon positions. For spherical nuclei, the collision ge¬ 
ometry is entirely a function of the impact parameter. 
However, in prolate deformed uranium, the geometry of 
the initial temperature distribution also depends on the 
relative spatial orientation of the two nuclei which can be 
described by the Euler angles between their long major 
axis. 

We focus in this paper on two limiting cases for fully 
overlapping uranium collisions. At one extreme are the 
“tip-tip” events, defined to occur when the major axes 
of the nuclei lie parallel with the beam direction. Re¬ 
siding at the opposite limit are the “body-body” events, 
where the major axis of each nucleus is perpendicular to 
the beam direction. We are interested in answering the 
question of how and with what precision we can distin¬ 
guish experimentally between these configurations. Their 


conceptual importance is explained in [7]. 

A preliminary account of parts of this work can be 
found in [13]; however, all results shown in Sec. IIIB 
involving hydro dynamically evolved events with ZDC se¬ 
lection constraints are new, and some of the results pre¬ 
sented in Sec. Ill A are based on larger hydrodynamic 
event samples. 


II. THE MODEL 

To model the initial energy density distribution of 
U+U collisions we employ the two-component (wounded 
nucleon/binary collisions) Monte-Carlo Glauber model. 
We use the deformed Woods-Saxon [14] distribution 

P(u 0 , +) = - - e ( r _ r (0 j¥ ,))/ d l 1 ) 

to sample the positions of nucleons inside a uranium nu¬ 
cleus. In Eq. (1) we take for the surface diffusiveness 
parameter d = 0.44 fm and for the saturation density 
parameter p 0 = 0.1660 fm -3 [15, 16]. The spatial config¬ 
uration of a uranium nucleus is deformed; we model its 
charge radius as [17] 

oo l 

r(0,<p) = r 0 (l + Yi X PirnYr(9, (/))), (2) 

1=1 rn= — l 

where 7+ = 6.86 fm is the average nuclear charge radius 
[16]. We assume the uranium nucleus is azimuthally sym¬ 
metric and choose the non-vanishing deformation param¬ 
eters (320 = 0.28 and /3 40 = 0.093 for the quadrupole and 
hexadecupole deformations along its main axis [18]. The 
choices of these parameters agree well with a recent re¬ 
analysis in [19], except for ^20 for which [19] gives the 
value 0.265. 

We use the Woods-Saxon density (1) to Monte-Carlo 
sample the nucleon centers and represent each nucleon in 
the transverse plane by a gaussian areal density distribu- 
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tion about its center: 

^=2^ e_ri/(2S) - (3) 

The width parameter B = ^nn(\/ s nn)/ 14.30 depends on 
collision energy as described in [ 20 ]. 

The sum of these gaussian nucleon density distribu¬ 
tions represents the nuclear density distribution for the 
sampled nucleus at the time of impact and is used to 
compute the initial energy density distribution gener¬ 
ated in the collision. For this calculation, we use the 
two-component Monte-Carlo Glauber model which adds 
contributions from binary collisions Nb and wounded nu¬ 
cleon participants N p with an adjustable relative weight 
[21]. The collision criterium that identifies binary colli¬ 
sions and wounded nucleons is evaluated probabilistically 
with the Gaussian nucleon profile (3) [20]. 

The binary collision term counts the entropy deposited 
by pairs of colliding nucleons and is modeled by a gaus¬ 
sian distribution with the same size as a nucleon (see 
Eq. (3) [22]; the total binary collision density per unit 
transverse area is 

n B c(r ± ) = Y^li,j^^e-\ F± - R ^\ /(2B) (4) 

h3 

where the sum is over all pairs of colliding nucleons. The 
normalization 7 \ is a T-distributed random variable with 
unit mean that accounts for multiplicity fluctuations in 
individual nucleon-nucleon collisions. 

Each struck nucleon is said to be wounded by (or 
participating in) the collision and contributes a portion 
of the initial entropy density distributed symmetrically 
about its center; the resulting total wounded nucleon 
density per unit area is given by 

n W N{ ? j _) = ^ e_|fi ‘ fi,i|!/(2B) ( 5 ) 

i 

where the sum is over all wounded nucleons in both nuclei 
and 7 i is again a fluctuating factor with unit mean. 

As mentioned, we model multiplicity fluctuations in 
single nucleon-nucleon collisions using 7 ij and 7 which 
are taken to be T-distributed random variables with unit 
mean and and with variances controlled by parameters 
Obc and Own, respectively. The generic T distribution 
with unit mean and scale parameter 6 is given by: 

1 / 0-1 -7/0 

r(7; ^ )= r(l/fl)6>v» ’ T^I 0 ’ 00 ) (6) 

The multiplicity fluctuations from wounded nucleons and 
binary collisions are related by requiring [ 22 ]: 

1 — a 

Qpp — 2 —® WN = a ® BC ' (?) 

The parameter 0 PP = 0.9175 was fit to multiplicity dis¬ 
tributions measured in p+p collisions, assuming multipli¬ 
city fluctuations coming purely from the initial state [23] . 


The distribution in the transverse plane of the de¬ 
posited entropy per unit volume is determined by mixing 
the binary collision and wounded nucleon sources using 


so(rj_) = 


TO 



a 


nwN(?±) + 


( 8 ) 


where To is the starting time for the (hydro)dynamical 
evolution of the collision fireball. We choose k s = 17.16 
and the mixing ratio a = 0.12 to reproduce the mea¬ 
sured charged multiplicities and their dependence on 
collision centrality in Au+Au collisions at 200 AGeV. 
The shape of the resulting energy density distribution 
in the transverse plane is calculated from the entropy 
density using the equation of state (EoS) s95p-vO-PCE 
from Lattice QCD [24]. The initial energy profile is 
evolved using the viscous relativistic fluid dynamic code 
package iEBE-VISHNU [22] with specific shear viscosity 
77 /s = 0.08. Simulations begin at time To = 0.6 fm/c and 
decouple at a temperature T^ ec = 120 MeV. The single 
particle momentum distribution is then computed using 
the Cooper-Fyre Formula. A full calculation of charged 
hadron observables that includes all hadronic resonance 
decay processes on an event-by-event basis is numeri¬ 
cally costly; for this reason we computed only the di¬ 
rectly emitted positively charged “thermal pions”, 7 r + 
and take this quantity as a measure for total charged 
multiplicity. At a fixed freeze-out temperature of 120 
MeV, the two quantities are related by a constant factor: 
dN ch /dr] ~ 4.6 dN n +/dy. 

The initial energy density profiles fluctuate from event 
to event. Each profile can be characterized by the r n - 
weighted eccentricity coefficients £ n and their associated 
“participant plane angles” 4> n : 


; in*„ _ f d? ± r n e in *e(? ± ) 
f dr±r n e(r±) 


(n > 2). (9) 


where (r, ip) are the standard polar coordinates in the 
transverse plane and e(?j_) is the initial energy density 
[25]. Through the hydrodynamic evolution, these spa¬ 
tial eccentricities {e n ,4> n } translate themselves into the 
anisotropic flow coefficients {u n ,T n } [26-28]: 


v , = v e in* n = f PTdp T d(f p e in ^dN/(prdprdipp) 
f PTdp T dippdN/ (p T dp T d(f p ) 

( 10 ) 

Apart from the Monte-Carlo Glauber model, there ex¬ 
ist various other initialization models. These include the 
IP-Glasma model [29], the MC-KLN model [30-33], and 
the TRENTO model [34]. As we will see, U+U collisions 
can provide experimental measurements to distinguish 
between these various initializations. 
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Figure 1. Panels (a,b) show the event-averaged eccentricities e n , before hydrodynamic evolution, panels (c,d) the event-averaged 
flows v n after hydrodynamic evolution. The left panels (a,c) represent 35,000 minimum bias events that include multiplicity 
fluctuations whereas the right panels (b,d) were obtained from a different set of 35,000 multiplicity-selected events, generated 
earlier without including multiplicity fluctuations, covering only the 0-5% centrality range. 


III. CONSTRAINING COLLISION GEOMETRY 
WITH MULTIPLICITY, FLOW, AND ZDC CUTS 

A. Eccentricity and flow coefficients as a function 
of multiplicity 

In Fig. 1, we present the centrality dependence of the 
initial eccentricities and the final anisotropic flow coef¬ 
ficients of thermal pions for harmonic order n = 2 — 5 
in U+U collisions at 193 AGeV. In Figs. la,c minimum 
bias results are shown as functions of the thermal pion 
yields, dN n + /dy. We notice that the variance of 52,4 and 
t > 2j 4 in “most central” (i.e. highest multiplicity) collisions 
are larger than in the rest of the centrality range. This 
is because in the most central collisions the two uranium 
nuclei are colliding with impact parameter b ~ 0 but, as 
a result of the large spatial deformation, not always with 
full overlap. A mixture of tip-tip and body-body col¬ 
lisions in these highest multiplicity events increases the 
variance of the initial £ 2,4 which in turn drives a larger 
variance in ^ 2 , 4 . 

In Figs. lb,d we focus on the 0-5% most central U+U 


collisions and increase the statistics to 35,000 events for 
just this bin. We find a “knee” structure in the high 
multiplicity regime (< 0.5% centrality) for both £ 2 and 
V 2 • This can be understood as follows: First, while the 
ellipticity in the transverse plane for a tip-tip collision 
is small (as the overlap area is approximately circular), 
body-body collisions produce ellipsoidally deformed ove- 
lap regions with larger ellipticities £ 2 . Second, although 
fully overlapping tip-tip and body-body collisions share 
the same number of participants, more binary collisions 
between nucleons can happen in the optically thicker 
tip-tip event, implying (in our two-component Glauber 
model) a larger initial dS/dy deposited for the tip-tip 
configuration. In the presence of fluctuations which lead 
to a range of £2 values for a given dS/dy and vice-versa, 
the larger average multiplicity in tip-tip collisions implies 
an increasing bias toward small £2 when selecting events 
with larger and larger values of dS/dy. This preferen¬ 
tial selection of tip-tip orientations at high multiplicities 
accounts for the appearance of a knee structure in the 
initial ellipticity [11] (Fig. 1). We see in Fig. lc that the 
knee is preserved after an event-by-event hydrodynamic 
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Figure 2. The ellipticity €2 as a function of dS/dy from the 
MC-Glauber model, for collisions roughly in the 0-5% cen¬ 
trality range, with (blue dashed line) and without (black solid 
line) multiplicity fluctuations from single p+p collisions. 


simulation when plotting the elliptic flow of the final par¬ 
ticle distribution as a function of multiplicity. 

We emphasize that that experimental results from 
STAR do not show this knee structure [12]. Consider¬ 
ing the preservation of the structure after hydrodynamic 
evolution as seen in Fig. 1, we conclude that, in contrast 
to Au+Au collisions where it has been extensively tested, 
the two component MC-Glauber model fails to correctly 
identify entropy production in ultra-central U+U colli¬ 
sions where the knee is predicted by the model but not 
found experimentally. Hence the non-linear dependence 
of multiplicity on the number of wounded nucleons ob¬ 
served in spherical Au+Au and Pb+Pb collisions as a 
function of collision centrality cannot be attributed to a 
binary collision component as implemented in the two- 
component MC-Glauber model. 

The IP-Glasma model [35] and the wounded- 
constituent-quark Glauber model [36] provide alternative 
explanations of the measured centrality dependence of 
the charged multiplicity without invoking a binary col¬ 
lision component. Gluon saturation physics as imple¬ 
mented in the IP-Glasma model is able to simultaneously 
accomodate a strong nonlinearity of dN/dy as a function 
of 7Vp art in Au+Au and Pb+Pb collisions and a weak 
dependence of dNjdy on collision orientation in central 
U+U collisions with a fixed number of participants, while 
the MC-Glauber model cannot [35]. Alternatively, it has 
been suggested in [36] that a model that produces en¬ 
tropy according to the number of wounded valence glu¬ 
ons (rather than wounded nucleons) can also reproduce 
the observed nonlinearity of dN/dy as a function of par¬ 
ticipant nucleons in Au+Au and Pb+Pb at RHIC and 
LHC, without a binary collision component. It would be 
interesting to study the prediction of such a model for 
central U+U collisions of varying orientations. 


Some corrections to the entropy production in these 
ultra central events arise from the inclusion of p+p mul¬ 
tiplicity fluctuations. In Fig. 2 we show that accounting 
for multiplicity fluctuations increases the average ellip¬ 
ticity of the collision fireball and weakens but does not 
erase the knee structure in £2 vs. dS/dy. Hence, this ef¬ 
fect alone does not appear sufficient to reach agreement of 
the MC-Glauber model with data for ultra-central U+U 
collisions. We acknowledge that more drastic fluctuation 
models [37] have been suggested in order to more success¬ 
fully adjust the theoretical predictions of MC-Glauber to 
experimental results. 


B. Selecting high overlap events with combined 
ZDC and multiplicity cuts 

In an experimental analysis of relativistic heavy ion 
collisions, the charged hadron multiplicity, dN^/dy and 
the elliptic flow coefficient V 2 can be used to classify 
events. Hydrodynamic studies have shown that the ini¬ 
tial £2 maps linearly to the V 2 of hadrons [28] and that 
the initial dS/dy is monotonically related to the final to¬ 
tal particle multiplicity, dN/dy [22]. Hence one should 
be able to use dS/dy and £2 from the initial conditions 
as proxies for the (computationally intensive) charged 
hadron dN^/dy and V 2 when testing our ability to se¬ 
lect the fully overlapping tip-tip and body-body U+U 
collisions. In this subsection we test this idea. 

In our analysis we make theoretical approximations for 
the use of experimental forward and backward zero de¬ 
gree calorimeters (ZDCs). Placed at zero degrees far 
from the colliding pair, ZDCs catch information about 
the spectator neutrons that pass through a collision with¬ 
out getting wounded. We classify our collisions by using 
the number of spectators N s = 476 — A part (where A part 
is computed from the MC Glauber model) to mimic the 
experimental ZDC signal [8]. For this study we look at 
65,000 intial condition events in the 1% most partici¬ 
pating ZDC range (N s < 19) and compare statements 
made on the basis of their ellipticities £2 and initially pro¬ 
duced entropy dS/dy with analogous statements made 
on the basis of elliptic flow V 2 and final charged multi¬ 
plicity dNch/dy for a subset of 25,000 events that were 
evolved using viscous relativistic fluid dynamics. The cal¬ 
culations include T-distributed multiplicity fluctuations 
in nucleon-nucleon collisions. Selecting the most partic¬ 
ipating ZDC collisions allows for a restriction of the set 
of collisions to more fully overlapping events. In such a 
regime, any initial geometric effects should come almost 
exclusively from the deformed shape of the uranium nu¬ 
cleus. 

We define the tip-tip and body-body event classes us¬ 
ing the pair of angles (# 1 , 2 ? ^ 1 , 2 ) from the two incoming 
nuclei, where 6 is the polar angle between the long major 
axis of the uranium nucleus and the beam direction and </> 
is its azimuthal angle in the transverse plane. An event 
is classified “tip-tip” if yj (cos 2 0 \ + cos 2 62 )/2 > a/3/2, 
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and “body-body” if both yj (cos 2 6 \ + cos 2 0 2 ) /2 < 1/2 
and |0i — 0 2 I < 7r/10. The polar angle constraints are 
chosen such that if the first nucleus has 6 1 = 0 the colli¬ 
sion is called “tip-tip” for 62 < 7r/4. Accordingly, if the 
first nucleus has 0\ = 7r/2, the collision is called “body- 
body” for 62 > 7r/4. These polar angle constraints imply 
that if 6 \ = O 2 , the common angle is less than 7 r/6 for 
“tip-tip” and greater than 7 t/ 3 for “body-body”. The 
additional azimuthal constraint included for body-body 
events is intended to force good alignment of the long 
major axes. 



Figure 3. Probability distributions for dS/dy (top) and 
dN^+fdy (middle, scaled by (dS / dy) / (dN 7T + / dy)) showing 
tip-tip and body-body collision contributions, within a se¬ 
lection of 1% ZDC events. The bottom panel shows the rela¬ 
tive probabilities for tip-tip and body-body events among all 
events of a given dS/dy. 


In Fig. 3 we plot the probability distributions for dS/dy 
and dN^+fdy. Using our collision definitions we can di¬ 
rectly read from the figure the likelihood of selecting a 
certain orientation based on a given multiplcity cut. Not¬ 
ing the evident congruence of the plots for dS/dy (upper) 
and dN n +/dy (middle), we rely on dS/dy in the bot¬ 
tom plot because of the improved statistics. We see that 
by cutting (within our 1% ZDC sample) on events with 
large dS/dy we can enrich the fraction of tip-tip events to 
about 50%, whereas cutting on low dS/dy enriches the 
fraction of body-body events, but never to more than 
about 20%. The 20% limit arises from admixtures from 
imperfectly aligned collisions that are not really ”full 
overlap”. The enrichment of tip-tip or body-body by 
varying dS/dy relies on the assumed two-component na¬ 
ture of entropy production which also produced the knee 
structure discussed before. Indeed, selection efficiency of 
specific collision geometries by cutting on dS/dy is model 
dependent. 

We therefore consider “event engineering”, an ap- 
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Figure 4. Probability distributions for £2 (a), q 2 (b) and V 2 (c) 
for different event classes within a sample of 1% ZDC events. 
The top half of each panel panels shows the distributions of 
tip-tip and body-body collisions, while the bottom half shows 
the relative contributions from each class. 


proach for selecting events by the magnitudes ^2 and V 2 of 
their Q 2 or V 2 flo w vectors. For high multiplcity events, 
Qi « VdNch/dy v 2 [38]. In Fig. 4 we plot the probabil¬ 
ity distributions for £2 (a) and compare them with those 
for q 2 (b) and V 2 (c). Both ^2 and V 2 are computed for 
directly emitted pions only, assuming that ~ Uf\ 
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and similar for q 2 . Although not quite true, this ap¬ 
proximation should cause at most a slight rescaling of 
the horiziontal axes in Figs. 4b, c [39], without changing 
the shape of the distributions. Calculating the full reso¬ 
nance decay chain for all 25,000 events would have been 
prohibitively expensive. Comparison of panels b and c 
shows very similar probability distributions for q 2 and 
v 2 , i.e. very little influence of fluctuations of the charged 
particle multiplicity dN^/dy on their shape. 

Since tip-tip events have on average smaller ellipticities 
(see upper panel in Fig. 4a), selecting events with small 
ellipticity (or, in experiment, small q 2 or v 2 ) enriches the 
tip-tip fraction. However, in this way we will never reach 
more than about 25% purity of the tip-tip sample. On 
the other hand, cutting the 1% ZDC events on large q 2 
will enrich the sample in body-body events, with a purity 
that can reach about 40% for the largest q 2 values. 

The current ZDC cut strategy can be further refined 
to increase the probability of selecting tip-tip events. 
Rather than looking at the ZDC signal in one of the two 
ZDC detectors or the sum of the ZDC signals in both 
detectors, we can look at the correlation of these two sig¬ 
nals. Events with roughly equal forward and backward 
ZDC signals (i.e. approximately equal numbers of spec¬ 
tators from both nuclei) provide a better definition of 
the categories full overlap, tip-tip, and body-body than 
events with asymmetric ZDC signals where all spectators 
come from only one of the colliding nuclei. The differ¬ 
ence in participants A7V part = | Apart, l - N part , 2 | quanti¬ 
fies the ZDC correlation in our model. Low values of A N p 
correspond to the most correlated forward and backward 
ZDC signals. To demonstrate one application, we recon¬ 
sider the purity fractions shown in the bottom part of 
Fig. 4b and now select from the sample only events in 
the lowest 25% of AA part . The selection on small values 
of A7V part eliminates from the sample asymmetric con¬ 
figurations that we loosely describe as “tip-body”. Col¬ 
lisions of this type produce low values of e 2 (hence q 2 ) 
without the angular criteria necessary to be considered 
tip-tip and therefore dilute the contribution of the true 
tip-tip configurations in the lower range of £2 (# 2 )- We 
show in Fig. 5 that selecting the lowest 25% of A7V part 
increases the selection efficiency of q 2 for tip-tip configu¬ 
rations by a factor of about 1.4. As a final comment, we 
point out that it might also be interesting to use ZDC 
correlations in the opposite way so as to select and study 
events with asymmetric tip-body configurations. 


IV. CONCLUSIONS 

Within the two-component MC-Glauber model for ini¬ 
tial energy production, the prolate deformation of the 
uranium nucleus was shown to generate a knee in the 
centrality dependence of the ellipticity of the initial tem¬ 
perature distribution. The knee was seen to be preserved 
by hydrodynamic evolution, after which it manifests itself 
in the centrality dependence of v 2 . Such a knee struc- 



Figure 5. The solid curve shows the distribution of events 
scaled according to their contribution to the total probability 
distribution for q 2 as seen in the bottom right panel in Fig. 4. 
The dashed curve shows the increased contribution of tip-tip 
collosions within the 25% of events having the smallest differ¬ 
ence in participants AV part (a proxy for ZDC correlation). 


ture is not seen in the STAR data. This rules out the 
two-component MC-Glauber model for initial energy and 
entropy production. An enrichment of tip-tip configura¬ 
tions by triggering only on high-multiplicity in the U+U 
collisions thus does not work. 

To increase the selection capability between different 
collision geometries, we impose combined cuts on initial 
conditions using the spectators (ZDC), dN/dy , and v 2 
(or q 2 ). For 1% ZDC events, we found that we could 
enrich tip-tip collision geometries to about 50% by cut¬ 
ting on high multiplicity within that sample, and body- 
body configurations to about 20% purity by selecting 
low-multiplicity events. These numbers include effects 
from multiplicity fluctuations, but they rely on the binary 
collision admixture in the two-component MC-Glauber 
model and are thus mo del-dependent. Since the two- 
component Glauber model is experimentally disfavored 
by the failure to observe the “knee” structure predicted 
by that model in a plot of v 2 vs. dN^/dy, these purity 
factors may not be reliable. 

To eliminate the model-dependence just mentioned we 
also studied the efficiency of selecting different collision 
geometries by “event engineering”, i.e. by cutting on 
q 2 or v 2 . In this case events selected for high q 2 can 
enrich body-body collisions to about 40% purity while 
cutting on low q 2 gives a tip-tip sample with about 25% 
purity. The latter can be boosted to about 35% purity by 
eliminating events with asymmetric ZDC signals. These 
results should not be sensitive to the binary collision ad¬ 
mixture in the two-component MC-Glauber model and 
thus should be less model dependent. 
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